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Abstract 

Here we present a method to find elementary first integrals of rational second 
order ordinary differential equations (SOODEs) based on a Darboux type procedure 
[16, 12, 13]. Apart from practical computational considerations, the method will 
be capable of telling us (up to a certain polynomial degree) if the SOODE has an 
elementary first integral and, in positive case, finds it via quadratures. 
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1 Introduction 



The differential equations (DEs) are the most widespread way to formulate the 
evolution of any given system in many scientific areas. Therefore, for the last three 
centuries, much effort has been made in trying to solve them. 

Broadly speaking, we may divide the approaches to solving ODEs in the ones 
that classify the ODE and the ones that do not (classificatory and non-classificatory 
methods) . Up to the end of the nineteenth century, we only had many (unconnected) 
classificatory methods to try to deal with the solving of ODEs. Sophus Lie then 
introduced his method [1, 2, 3] that was meant to be general and try to solve any 
ODE, i.e., non-classificatory. Despite this appeal, the Lie approach had a short- 
coming: namely, in order to deal with the ODE, one has to know the symmetries 
of the given ODE. Unfortunately, this part of the procedure was not algorithmic 
(mind you that the classificatory approach is algorithmic by nature). So, for many 
decades, the Lie method was not put to much "practical" use since to "guess" the 
symmetries was considered to be as hard as guessing the solution to the ODE itself. 
In [4, 5], an attempt was made to make this searching for the symmetries to the 
ODE practical and, consequently, make the Lie method more used. 

Even thought the attempt mentioned above was very successful, the procedures 
applied to find the symmetries were heuristic. So, a non- classificatory algorithmic 
approach was still missing. The first algorithmic approach applicable to solving 
first order ordinary differential equations (FOODEs) was made by M. Prelle and 
M. Singer [6]. The attractiveness of the PS method lies not only in its basis on 
a totally different theoretical point of view but, also in the fact that, if the given 
FOODE has a solution in terms of elementary functions, the method guarantees that 
this solution will be found (though, in principle it can admittedly take an infinite 
amount of time to do so). The original PS method was built around a system of 
two autonomous FOODEs of the form x = P(x,y), y = P(x,y) with P and V in 
C[x,y] or, equivalently, the form y' = R(x,y), with R(x,y) a rational function of 
its arguments. 

The PS approach has its limitations, for instance, it deals only with rational 
FOODEs. But, since it is so powerful in many respects, it has generated many 
extensions [7, 8, 9, 10, 11, 12, 13, 14] 

Nevertheless, all these extensions deal only with FOODEs. In particular, the 
second order ordinary differential equations (SOODEs) play a very important role, 
for instance, in the physical sciences. So, with this in mind, we have produced [17] 
a PS-type approach to deal with SOODEs. This approach dealt with SOODEs that 
presented elementary 2 solutions (with two elementary first order invariants). 

Here, we present a different approach that, besides dealing with a much broader 
class of SOODEs (those with at least one elementary first order invariant), does not 
depend on a conjecture about the general structure of the first order invariants. 

In section 2, we present the state of the art up to the present paper. In the 
following section, we introduce some important theoretical results for the building 
of the algorithm to find the integrating factor. In section 4, we present the algorithm 
for finding the integrating factor with examples of its application. Finally, we present 
our conclusions and point out some directions to further our work. 

2 For a formal definition of elementary function, see [15]. 
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2 Earlier Results 

In the paper [6], one can find an important result that, translated to the case of 
SOODEs of the form 

// M(x,y,y') , 

y = ^77 k = nx, y, y ), (i) 

N(x,y,y') 

where M and N are polynomials in (x,y,y'), can be stated as: 

Theorem 1: // the SOODE (1) has a first order invariant that can be written 
in terms of elementary functions, then it has one of the form: 

m 

I = wo + ^Qln(wj), (2) 

i 

where m is an integer and the w's are algebraic functions 3 of (x,y,y'). 
The integrating factor for a SOODE of the form (1) is defined by: 

R(^-y") = d -^^ (3) 

where ^ represents the total derivative with respect to x. 

Bellow we will present some results and definitions (previously presented on [17]) 
that we will need. First let us remember that, on the solutions, dl = I x dx + I y dy + 
I y i dy' = 0. So, from equation (3), we have: 

R((f> dx - dy') = I x dx + I y dy + I y < dy' = dl = 0. (4) 
Since y' dx = dy, we have 

R[{4> + Sy')dx-Sdy-dy']=dI = 0, (5) 

adding the null term S y' dx — S dy, where S is a function of (x,y,y'). From equation 
(5), we have: 

I x = R(4> + Sy'), 

I y = -RS, (6) 

Iy' = —R, 

that must satisfy the compatibility conditions. Thus, defining the differential oper- 
ator D: 

D = d x + y'd y + (t>d yl , (7) 
after a little algebra, that can be shown to be equivalent to: 

D[R] = -R(S + <f> yf ), (8) 
D[RS] = -R4> y . (9) 



3 For a formal definition of algebraic function, see [15]. 
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3 New theoretical results concerning the func- 
tion S 



Let us start this section by stating a corollary to theorem 1 concerning S and R. 

Corollary 1: If a SOODE of the form (1) has a first order elementary invariant 
then the integrating factor R for such an SOODE and the function S defined in the 
previous section can be written as algebraic functions of (x,y,y'). 



Proof: Using the above mentioned result by Prelle and Singer, there is always 
a first order invariant I = wq + YT c « ^ n ( w i) f° r the SOODE. So we have, using 
equation (3), 

- y") = h + y'ly + y"i y > ^R = -iy (10) 

where I u = d u I. From equation (2), we have: 



m 

Wi,.i 



LU i 

i 1 

Then I y i is an algebraic function of (x,y,y') and, by equation (10), so is R. 
From equations (6), one can see that: 



Therefore, S is also an algebraic function of (x,y,y').D 



Besides that, working on equations (8) and (9), we get [17]: 

M 

AT 



D[S] = S 2 + ( f>y'S-<i )y = ^, (13) 



where M. and N are given by 

M = {NSf + (NMyi - MNyi) S - (NM y - MN y ), (14) 
M = N 2 . (15) 

Concerning eq.(13) we can demonstrate the following theorem: 

Theorem 2: Consider the operator defined by D$ = M ds +N D. If P is an 
eigenpolynomial of Ds (i.e., Ds[P] = XP, where X is a polynomial) that contains S, 
then P = defines a particular solution of eq.(13). Conversely, If P is a polynomial 
that contains S, such that P = defines a particular solution of eq.(13), then 
P is either an eigenpolynomial of D$ or P is an absolute invariant of the Lie 
transformation group defined by Dg- 



Proof: In order to demonstrate theorem 2 we will, first, prove the following lema: 
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Lema 1: If eq.(13) (D[S] = M./M) has an algebraic solution defined byJ2i a iS % = 
0, where the ai are polynomials in (x,y,y'), then 

Y (AfDfa]? + MaiiS^ 1 ) = 0. (16) 

i 

Conversely, if J2i(N 'D[ai]S l + MaiiS' l ~ l ) = then J2i a iS l = defines an 
algebraic function (S) as a particular solution of eq.(13). 

Proof of Lema 1: We begin by proving the first part of the Lema. Let 
J2i a iS l = define an algebraic particular solution of eq.(13). Applying the op- 
erator D on J2i a iS % = 0) one g e t s: 

Y (d^S* + aiiS^DlS]) = o Y (^M^ + w^jp) = °- 

i i 

Multiplying this by J\f we get eq.(16). 

Let us now prove the converse. Consider that eq.(16) applies. Multiplying it by 
N and remembering that S obeys D[S] = M/Af, one gets: 

Y (d^S* + OiiS*- 1 ^) = =► - £ (D[a t ]S*) /£ (oiiST 1 ) = ^ (17) 

i i i 

However, note that, if we have an algebraic function (B) defined by J2i hB l = 0, 
then D[B] can be found as follows: 



D 



^Y ( D M Bi + &ii-B i_1 £)[B]) = 



D[B] = -Y(D[b i ]B i )/Y( b i iBi ~ 1 ) 



Therefore, one can see that eq.(17) can be put in the form: 



D[S] = (18) 



where S is an algebraic function defined by J2i a i^ % = 0. □ 



Now, using Lema 1, we are going to demonstrate Theorem 2. 

Consider that P is an eigenpolynomial of D$ that contains S. By definition, 
Ds[P] = XP, where A is a polynomial. So, writing P = J2i a iS l , where the Oj's are 
polynomials in (x,y,y'), we have: 

Y^Dla^S 1 + Ma i iS i - 1 ) = XP (19) 

i 

Over the algebraic function defined by P = 0, we get J2i {■N'D[ai]S % + MaiiS 1 ^ 1 ) 
0. This, using Lema 1, implies that J2i a iS l = defines a particular solution to 
eq.(13). 
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Conversely, consider that P is a polynomial that contains S, such that P = 
J2i a iS l = defines a particular solution of eq.(13). Again, via Lema 1, we have 
that: 

(AfDladS* +Ma i iS i - 1 ) = 0. (20) 

i 

Consider now the following operator: 

= N\x 1 ,x 2 ,x 3 ) (d Xl +x 3 d X2 + ^ Xl ' X2,X3 ] d X3 ) +M(x 1 ,x 2 ,x 3 ,x A )d XA (21) 
V N(x 1 ,x 2 ,x 3 ) ) 

Applying O to a generic polynomial V = J2i Cix\, where the Cj's are polynomials in 
(xl,x2,xS), one obtains: 

0[V] = J2(N 2 (x 1 ,x 2 ,x 3 )V[c i ]xi + M(x 1 ,X2,X 3 ,X4)Ciix\ *) (22) 

i 

where V = (d xi + x 3 d X2 + %$^d X3 ). 

Since the terms multiplying the partial derivatives are polynomials, applying O 
to a polynomial will generate a polynomial. So, from eq.(22), we have: 

[N 2 (xi,X2,x 3 )V[ci]x\ + M(x 1 ,x 2 ,x 3 ,X4)c i ix l 4 ~ 1 ^ = Q (23) 

i 

where Q is a polynomial in (x±, x 2 , x 3 , X4). 

Note that the left-hand side of equations (23) and (20) are formally equiva- 
lent. Consider that the hypothesis of the theorem 2 apply, i.e., V = J2i c i x \ = 
defines a algebraic function x^(xi, x 2 , x 3 ) that is a particular solution of -D[x4] = 
Ai(xi, x 2 , x 3 , Xi)/Af(xi, x 2 , x 3 ). Lema 1 implies then that, if we substitute the func- 
tion X4 = RootOf(J2i Cix\) into eq.(23), we get = Q. Therefore, V = =>■ Q = 
0. Since both V and Q are polynomials we have two possibilities: 

• Q is identically null - in that case, V is an absolute invariant of the Lie trans- 
formation group defined by the operator O. 

• V is a factor of Q - in that case, Q = XP, where A is a polynomial in 
(x±, x 2 , x 3 , X4) and Q is a relative invariant of the Lie transformation group 
defined by the operator O. 

This completes the proof of theorem 2.D 



From the results above we may finally conclude: 

Corollary 2: If the SOODE (1) has an elementary first integral, then there 
is a polynomial P ( containing S ) that is either an eigenpolynomial of D$ or is an 
absolute invariant of the Lie transformation group defined by D$- 



Proof: From corollary 1, since eq.(l) has a first order elementary invariant, 
there is an algebraic function S that satisfies eq.(13). So, by definition, there is a 
polynomial P (that contains S) such that P = defines S. From theorem 2, this 
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implies that the polynomial P is either an eigenpolynomial of Ds or is an absolute 
invariant of the Lie transformation group defined by D$.0 

These theoretical results provide us an algorithm to find S. Briefly, in words, 
what we have to do is the following: 

• Find the eigenpolynomials of the Ds operator containing S. Let us call them 
Pi- 

• In order to find S, we have to choose one of those Pi (let us call it P) that 
contains S and solve the equation P = for S. 

• If you succeed in doing the above steps, you would have found S for the 
SOODE you are dealing with. 

The importance of these results and the above sketched method is that, as we 
shall briefly see, the finding of the algebraic function S will allow us to produce a 
semi-algorithmic method to find the integrating factor (and, consequently, the first 
order invariant) for SOODEs of the type described in eq.(l). 

4 Finding the integrating factor and the first 
integral 

In what follows, we are going to demonstrate a result concerning the general struc- 
ture of R, Let us state that result as a theorem. 

Theorem 3: Consider a SOODE of the form (1), where M and N are poly- 
nomials in (x,y,y'), that presents an elementary first order invariant I. Then the 
integrating factor R for this SOODE can be written as: 

R = Uf?> (24) 

i 

where fi are irreducible polynomials in (x,y,y',S), which are eigenpolynomials of 
the operator V = N D or are factors of N and are non-zero rational numbers. 

Proof: 

Suppose that the hypothesis of the theorem is satisfied. Re- writing equation (4): 

— [(M + S N y')dx - S N dy - N dy'] = dl = (25) 

For the sake of simplicity, let us establish some notation: M = (M + y' N S), 
N = -(N S) and R = (§). We can write I x = RM and I y = -RN and imposing 
the compatibility condition I xy = I yx , we have d y (RM) = —d x (RN). Expanding 
this 

R y M + R X N = -R(M y + N x ) = -(M y + N x ) (26) 

R 

where D = Nd x + Md y . 

Since M and TV are polynomials in (x, y, S) 4 , equation (26) is formally equivalent 
to the equation which establish the conditions for the theorem by Prelle and Singer 

Note that M and N are polynomials in y' as well. 
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[6] with the extension by Shtokhamer [7] 5 . More formally, consider the FOODE 
defined by: 

where M and N are the polynomials obtained by replacing y' by a constant k in 
M and iV respectively. Then we can write the solution for the FOODE (27) as 
I = C, where / is defined by replacing y' (by the constant k) on /. By hypothesis, 
/ is elementary and, consequently, / is also elementary. Therefore, by the Theorem 
due to Prelle and Singer [6] and the Shtokhamer's extension [7], the integrating 
factor R (defined by replacing y' (by the constant k) on R) for equation (27) can 
be written as Yli f,™* , where fi are irreducible polynomials in (x,y,S). We could 
proceed analogously for the other two possible pairings of variables (i.e., (x,y') and 
(y, if)). More explicitly, in the reasoning leading to equation (26), we could have 
imposed the compatibility condition I xy , = I y / X and, in the same fashion we did 
above for the pairing (x,y), conclude that the integrating factor can be written 
as a product of irreducible polynomials in (x,y',S). The same can be done for 
the pairing (y,y'). So, we can conclude that R can be written as a product of 
irreducible polynomials in (x,y,y',S). Since R = (#), and N is a polynomial in 
(x,y,y'), finally we have that R = FJi /"% where /, are irreducible polynomials in 
(x,y,y',S) and rii are non-zero rational numbers. 

Let us now prove that the fi are eigenpolynomials of the operator T> or factors 
of N. From equation (8) and remembering that R = we have: 

D\RN] N M y i — M N y i 

RN N 2 1 ' 

multiplying both sides of (28) by N 2 , we get 

N 2 ^1B- + ND[N] = -N 2 S- {N M yl - M N y >) (29) 
R 



and finally, 



AT 321 = -ND[N] -N 2 S- (NM y , - M N y ,). (30) 
R 



Since the right-hand side of (30) is a polynomial in (x, y, y', S), so is N m. Using 



(24) this implies that 



jv*B = 2>atH« (31 ) 



is a polynomial. Since fi are irreducible, independent polynomials, we may conclude 
that: either fi\V[fi] or fi is a factor of N.O 

In order to obtain the integrating factor R, we are going to use eq.(8). Dividing 
it by R, we get: 

omR--v + M --(s + ™£j™t ) (32) 



5 Notice that, since S is an algebraic function of (x,y,y'), it is a root of a polynomial equation of 
the form ^iPiS 1 = 0, where pi are polynomials in (x,y,y'). Therefore, its derivatives can be expressed 
in terms of rational powers of itself. So, one may consider that would be dealing with a Shtokhamer 
extension with U = S. 
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Multiplying both sides of the equation above by N 2 , one has: 



N 2 D[R]/R = -(SN 2 + NM y , - MN y >) =► NV[R]/R = -(SN 2 + NM y , - MN y ,) 

(33) 

Due to the results of theorem 3, the equation (33) above it can be solved in the 
same manner as in the methods inspired on the Prelle- Singer Procedure. Once R 
is found as a function of (x, y, y', S), we can substitute the known S (see section 3) 
and, finally, from eqs.(6), find the first order invariant via quadratures: 



I(x, y, y') = Jr (</> + Sy') dx - J '[RS + |- J R (<f> + Sy') dx]dy 



dy 

I R + ^ (J R (0 + s vl dx - J[ RS + J + s v') dx \ d y 



dy'. (34) 



That concludes the presentation of our proposed approach which is a semi- 
algorithm procedure to reduce soodes of the form of eq.(l). 

4.1 Example 

Here we are going to briefly introduce an example of a non- trivial 6 SOODE that 
illustrates the results displayed above. 
Consider the SOODE: 

v" = - (yf+lMl^ll (3 5) 

For this SOODE, the operator D becomes: 

D = {y' 2 + 2y'Z-l-2y')d y , + (x + y'x + y'y+l-y' 2 )y'd y + 

(x + y' x + y'y + l-y 12 ) d x , (36) 

and the operator D$ is: 

D S = ND + Md s (37) 

where M = 2 s 2 xy'y + 2 s 2 y' 2 xy + 7 sy' 2 x + 4 sy' 3 x + s 2 + s 2 x 2 + 2 s 2 x - 2 s 2 y' 2 + 
s 2 y' A - sx + 4 sy' 2 - 2 sy' 4 + sy + 4 sy' 3 y - 2 s 2 y' 2 x + s 2 y' 2 y 2 - 2 s 2 y' 3 y + 2 s 2 y'x 2 + 
2 sy'x - 2 s - y' + s 2 y /2 x 2 - 2 s 2 y /3 x + 2 s 2 y'y - 2 y' 2 + y' 3 + 2y' A + sy' 2 y + 2 s 2 y'x. 

The polynomial P = (x + y'x + i/'y + 1 — y' 2 ^j S + 1 — y' 2 is a eigenpolynomial 
of the -Ds operator. So, P = defines 5 as: 

S = - 2 (38) 

x + y'x + y'y + 1 - y' 

In order to calculate the integrating factor R, we need (see equations 31 e 32) 
to find the eigenpolynomials of the D operator. These are found to be (up to the 
first degree): 

h = 1+2/ 

h = i-y' 

h = l + 2y' 



;> For instance, the MAPLE package, release 7, could not reduce it. 
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and, for this particular case, these are sufficient (together with N) to build the 
integrating factor R. Combining equations 31 and 32 and solving for the raj's we 
have that {rai = — 3/2, n 2 = — 3/2, ra 3 = 0}. So, R is found to be 7 : 

r = * + y' x + y'y + l -£ m 

((y'-l)(l + y')) 3/2 
This leads to the first order invariant given by: 



I= X -p^ + lJy' + 4^?). (40) 

V-i + y' 2 



5 Computational considerations 

Here, we are going to introduce an alternative way of calculating the integrating 
factor for a given SOODE. Why do we call this section computational considera- 
tions? The reason is that the material contained in this section is not applicable 
to all SOODEs 8 of the form of eq.(l) and, indeed, we do not have a criteria of 
applicability But, on the other hand, the method we are about to present is "less 
costly" computationally, therefore could be useful in practical applications. In the 
following, we will be using some results from the Lie theory. 

Let us show a relation between the function S and a symmetry of the SOODE 
(1). Making the following transformation 



S— £M (41) 



eq.(13) becomes 



D 2 [r)} = <f>y,D[ri] + <t> y n. (42) 

Erom Lie theory we can see that eq.(42) represents the condition for a SOODE 
(1) to have a symmetry [0,7/]. So, from (41) we can find a symmetry given by 

rj = e -fSdx (43) 

Looking at eq. (5) we can infer that R is also an integrating factor for the auxiliary 
FOODE defined by 

where x is regarded as a parameter. Besides, for this FOODE, [77, D[rj\] is a point 
symmetry. So, R is given by 

If in (45) we use r\ defined by (43) we would get a singular R. However, eq.(42) 
has another solution independent from 7/ given by 



7 Note that R is formed by the product /f 1 / 3 " 3 iV 1 = f{ 3/2 f 2 3/2 N 1 . 

8 Note that the method presented on the previous section is a general semi-algorithmic approach to 
deal with SOODEs of the form (1) that present an elementary first order invariant. 
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r e J ,p v" tx 

rj = T) \ 7. — dx. (46) 

J rj z 

Using this rj in (45) we get R and, once this is done, we can calculate the first 
integral I by using simple quadratures. 

Actually, trying to solve eq.(43) could be complex. The operator J S dx that 
appears on equation (43) is meant to be the inverse of operator D (defined on 
(13)) 9 . As one might expect, actually finding rj from (43) could be trick. Actually, 
in general, it will be impossible to integrate (43) and that is the reason why this 
shortcut is not applicable in general. 



5.1 Example 

Consider the SOODE 

2y' 2 x 2 -2y' 2 -2xy'y-y 2 x A + 2y 2 x 2 -y 2 

y = M^T) • (47) 

Constructing the Ds operator we get that 

y 2 (x 2 -l) -S 2 y 2 -y' 2 + 2Sy'y (48) 

is an eigenpolynomial of it. Then S is given by S = y'/y + V x 2 — 1 and 77, rj are 
respectively 



y/x+y/x^T __\/x + Vx 2 -lf: 



'x—l\/x+le x V^ * 



v+Voc' 2 — 1 



dx 



ye 2 ye ~ 

R can be written as 



(49) 



y\ x + V x 2 — 1 , N 

R =— 7TT ( 5 °) 

VX — le 2 



leading to the following first integral: 



V x + Vx 2 - 1 (2y'y + Vx 2 - ly 2 ) 



r—, r ^x 2 -i 

Vr — le 2 



(51) 



6 Conclusion 

In [17], we have developed a method, based on a conjecture, to deal with SOODEs 
that presented an elementary solution (possessing two elementary first order invari- 
ants). 

In that same paper, we have introduced a function S to transform the Pfaffian 
equation related to the particular SOODE under consideration into a 1-form pro- 
portional to differential of the first order differential invariant. That function S was 
instrumental in finding the integrating factor for the SOODE. 

9 Actually, this meaning for the operator J is true for the whole of the present section. 
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Here, in the present paper, we introduce many theoretical results concerning 
that function S and present an way to calculated it via a Darboux-type procedure. 
Here also the function S is used to produce the integrating factor. 

We then demonstrate general results about the structure of the integrating factor 
R and that we can calculate it using a procedure very similar to the one inspired 
by the original work by Prelle-Singer [6] (applicable for FOODEs). These results 
assure that we have a semi- algorithmic procedure to deal with rational SOODEs, 
presenting at least one elementary first order invariant, i.e., there is no need to 
posses two such invariants and, consequently, we can cover a much broader class of 
SOODEs than before [17]. 

In the above section, we introduce an alternative way to calculate R from the 
knowledge of S for restricted cases. The motivation for that method is that, the 
general method, sometimes, can be computationally very demanding and it may be 
worth to have a go in the "shortcut" before embarking on long calculations (despite 
the general case being the sure think). 

We are searching for better algorithms for the many steps of the method pre- 
sented here in order to make it computationally more efficient. We are also working 
on a full computational implementation of the method as it stands today. 

References 

[1] Stephani, H. Differential equations: their solution using symmetries, ed. 
M.A.H. MacCallum, Cambridge University Press, New York and London 
(1989). 

[2] GW. Bluman and S. Kumei, Symmetries and Differential Equations, Applied 
Mathematical Sciences 81, Springer- Verlag, (1989). 

[3] P.J. Olver, Applications of Lie Groups to Differential Equations, Springer- 
Verlag, (1986). 

[4] E.S. Cheb-Terrab, L.G.S. Duarte and L.A.C.P. da Mota, Computer Al- 
gebra Solving of First Order ODEs Using Symmetry Methods. Com- 
put.Phys.Commun., 101, (1997), 254. 

[5] E.S. Cheb-Terrab, L.G.S. Duarte and L.A.C.P. da Mota, Computer Al- 
gebra Solving of Second Order ODEs Using Symmetry Methods. Corn- 
put. Phys.Commun., 108, (1998), 90. 

[6] M Prelle and M Singer, Elementary first integral of differential equations. 
Trans. Amer. Math. Soc, 279 215 (1983). 

[7] R. Shtokhamer, Solving first order differential equations using the Prelle-Singer 
algorithm, Technical report 88-09, Center for Mathematical Computation, Uni- 
versity of Delaware (1988). 

[8] C B Collins, Algebraic Invariants Curves of Polynomial Vector Fields in 
the Plane, Preprint. Canada: University of Waterloo (1993); C B Collins, 
Quadratic Vector Fields Possessing a Centre, Preprint. Canada: University of 
Waterloo (1993). 

[9] C. Christopher, Liouvillian first integrals of second order polynomial differential 
equations. Electron. J. Differential Equations 1999, No. 49, 7 pp. (electronic). 



12 



[10] C. Christopher and J. Llibre, Integrability via invariant algebraic curves for 
Planar polynomial differential systems Ann. Differential Equations 16 (2000), 
no. 1, 5-19. 

[11] J. Llibre, Integrability of polynomial differential systems, Handbook of Differ- 
ential equations, Ordinary Differential Equations, volume 1, Chapter 5, pages 
437-531. Edited by A. Canada, P. Drabek and A. Fonda. Elsevier B.V. 2004. 

[12] L.G.S. Duarte, S.E.S.Duarte and L.A.C.P. da Mota, A method to tackle first 
order ordinary differential equations with Liouvillian functions in the solution, 
in J. Phys. A: Math. Gen. 35 3899-3910 (2002). 

[13] L.G.S. Duarte, S.E.S.Duarte and L.A.C.P. da Mota, Analyzing the Structure 
of the Integrating Factors for First Order Ordinary Differential Equations with 
Liouvillian Functions in the Solution, J. Phys. A: Math. Gen. 35 1001-1006 
(2002) 

[14] L.G.S. Duarte, S.E.S.Duarte, L.A.C.P. da Mota and J.F.E. Skea, Extension 
of the Prelle-Singer Method and a MAPLE implementation, Computer Physics 
Communications, Holanda, v. 144, n. 1, p. 46-62, 2002 

[15] Davenport J.H., Siret Y. and Tournier E. Computer Algebra: Systems and 
Algorithms for Algebraic Computation. Academic Press, Great Britain (1993). 

[16] Y K Man and M A H MacCallum, A Rational Approach to the Prelle-Singer 
Algorithm. J. Symbolic Computation, 11 1-11 (1996), and refferences therein. 

[17] L G S Duarte, S E S Duarte, L A C P da Mota and J E F Skea, Solving second 
order ordinary differential equations by extending the Prelle-Singer method, J. 
Phys. A: Math.Gen., 34 3015-3024 (2001). 



13 



